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Abstract 

As well as known, the black hole evaporation problem is famous prob- 
lem. Because the S.W. Hawking found the black holes emit light at the 
future null infinity as a thermal radiation |T], we think that the black holes 
may be vanish. However, to prove this problem, we should solve field equa- 
tion, i.e. forth order partial differential equations [^[3]. However, we can 
find a method to solve this equation, and we could prove that the black 
holes finally vanish. To solve this problem we use dynamical horizons 
equation in the Vaidya spacetime. 



1 Introduction 

We can show the black holes finally vanish by the analytical method and by 
the dynamical horizon equations. As we proved paper at 2006 in |4], the black 
holes finally vanish, by solving the dynamical horizon equations, and present 
paper mainly based on the paper. Note that the dynamical horizon equations 
and the Einstein equations are same, and to solve dynamical horizon equations 
is same as to solve field equation. However, in the previous paper, we use the 
approximations that we use the negative energy [5] near the event horizon not for 
the dynamical horizons. So in this paper we do not use the approximations. Our 
paper explain the result of previous paper and we show black hole evaporation 
problem without any approximations. 

The black hole evaporation problem starts from S.W. Hawking [T] at 1975. By 
his study we could find the black holes emit lights as the thermal radiation. And 
we can think that the black holes may be vanished by this effect i.e. Hawking 
effect. However, the proof of this problem seems to be hard, because we should 
solve field equations, and the field equations are forth order partial differential 
equations with non-linear term. We can not find the mathematical methods to 
solve such equations. However, by the study of the dynamical horizons [6] [7], 
we can solve this problem. 
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We can explain shortly the previous works of the black hole evaporation. 
The orthodox work may be Page's result [8] as M oc — M~^, it comes from the 
dimensional analysis. The string approach is [5][Tn][II], or semiclassical theory 
typically using apparent horizon [T2]. Hiscok studied spherical model of the 
black hole evaporation using the Vaidya metric, which we also use in present 
work, to solve the black hole evaporation problem. However, he simply set a 
model not taking account of the field equation. Hajicek's work[14] treated the 
black hole mass more generally than our present case. However, he did not 
use the field equation either. One of the more recent studies is Sorkin and 
Piran's work [IS] on charged black holes. And neutral case has been done by 
Hamade and Stewart [TB]. Their conclusion is that black hole mass decreases or 
increases depending on initial condition. They used a model of the double null 
coordinates, and obtained a numerical result. But they did not consider the 
Hawking effect directly but they used massless scalar field as a matter. Bre- 
vik and Haines calculated primordial black hole evaporation [TJ. Very recently 
Hayward studied black hole evaporation and formation using the Vaidya metric 
[T5] . It seems no analytical equation has been proposed for the black hole mass 
with the Hawking effect taken into account. 

In a section [2] we explain our previous result [1] by short sentences. And 
in a section [3J we explain the calculations of the dynamical horizons equation 
without any approximations. Finally in a section HJ we discuss and conclude 
our results. 

2 Previous results by short sentences 

We can show our previous result in three subsections. In the subsection 12.11 
we explain the dynamical horizons equation. In the subsection 12.21 we explain 
the Vaidya spacetime. And in the subsection 12.31 we show our previous result 
by calculating the dynamical horizons equation. The subsection ?? and 12.31 is 
the review of the previous paper and mainly same as ??. The reason why we 
use the dynamical horizons equation is it only needs information of the horizon 
surface, and reason why we use the Vaidya spacetime is that the dynamics of 
the horizon is free by the fluent matter. 

2.1 Dynamical horizons 

The study of the Dynamical horizon is mainly depend on Ashtekar and Kr- 
ishnan's works [6] [7] , and we omit the main calculations to derive the dynamical 
horizons equation. The definition of the dynamical horizons are as follows. 

Definition. A smooth, three-dimensional, spacelike submanifold H in a 
space-time is said to be a dynamical horizon if it is foliated by preferred family 
of 2-spheres such that, on each leaf S, the expansion Gj;) of a null normal l"" 
vanishes and the expansion of the other null normal n" is strictly negative. 
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This definition is similar to the apparent horizons. For simplicity, the defini- 
tion of the black hole is that the black hole does not emit even light. Although 
the above definition was given by Ashtekar and Krishnan, the definition define 
only spacelike horizons. So we can modify the above definition, as follows. 

Definition (modified version). A smooth, three-dimensional, spacelike or 
timelike submanifold in a space-time is said to be a dynamical horizon if it 
is foliated by preferred family of 2-spheres such that, on each leaf S, the expan- 
sion of a null normal 1°" vanishes and the expansion 6(„) of the other null 
normal is strictly negative. 

The timelike dynamical horizon is also called timelike membranes. 

There is the dynamical horizons equation. From the Einstein equation on 
the dynamical horizons, we can use 3-1-1 decomposition and then 2+1 decom- 
position, and finally using the Gauss-Bonet theorem, we obtain 



Here, Ri is the radios of the dynamical horizons before the energy fiows in the 
horizons and i?2 is the radius of the dynamical horizons after the energy flows 
in the horizons. And G is the gravitational constant, and integration is on the 
dynamical horizon A.H , and Tab means stress-energy tensor with cosmological 
constant, and the f" is the normal timelike vector diagonal to horizons and 
is the Nfi — \dR\ where R is the black hole radius. And the second term of the 
integrand is the form of the Bondi energy. However, the second term become 
vanish because we only treat spherically symmetric spacetime i.e. the Vaidya 
spacetime. Then the equation ([Ij becomes as. 



In this paper we use this simple equation. We can easily calculate the timelike 
dynamical horizon equation only replacing f"" and f"^. And the way of 3-1-1 and 
2-1-1 decomposition is replaced. 

2.2 Vaidya spacetime 

The Vaidya metric is of the form 




(1) 




(2) 



Fdv^ + 2Gdvdr + r^dn^, 



(3) 
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where F and G are functions of v and r, and is null vector and r is the area 
radius, and M is the mass defined byM=|(l — a function of v and r. 
This metric is spherically symmetric. By substituting the Vaidya metric (2) 
into the Einstein equation so that we can identify the energy-momentum tensor 
Tab as 



2-^.- 

2G r 



rG 



(4) 
(5) 
(6) 



We do not need to check that the solution of the dynamical horizon equation sat- 
isfies the Einstein equation. Because we would like to consider the Schwarzschild 
like metric, we set v = t + r* , where r* is tortoise coordinate with dynamics 



r* = r + 2M (v) In 
For later convenience, we write. 



2M{v) 



1 



dr 
dr* 



(7) 



(8) 



There are two null vectors. 
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(9) 



corresponding to the null vector and the other is 
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(10) 



Here we multiply a "'^ so that 1°" = v"". This choice of the null vector 1°" is 
explained in figure [T] From now on we put. 



F = 
G = 



2M{v) 



1, 



(11) 
(12) 



in a similar form to the Schwarzschild metric, assuming that M{v) is a function 
of V only. For a constant M, the metric coincides with the Schwarzschild metric. 
We calculate the expansions 8(;) and 8(„) of the two null vectors because 
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Figure 1: For the case that the dynamical horizon decreases, we should choose 
= — + f" so that Z° points into the dynamical horizon. 



the definition of the dynamical horizon requires one of the null expansions to 
be zero and the other to be minus. The result is, 

= Ic^F-a) (13) 

ar\ -F + 2a J ^ ^ 

From 6(i) = we get, 

2F-a = 0. (15) 

we can check that the other null expansion 0(„)is strictly negative. There- 
fore in this case, we can apply the dynamical horizon equation. In the usual 

Schwarzschild metric with dynamics, both expansions become zero. This is the 
one of the reasons why we choose the Vaidya metric. By inserting equation (7) 
to equation (14), we obtain 



a = F (l - 2M.. m - l) + M(./2M-1) ^-) " 



(16) 



Note that a is proportional to F. Now we solve 0(;) = 0, to determine the 
dynamical horizon radius as 

2F-a = 2F-F( l-2M^ln( — -1 | + — -M^ | =0. (17) 

V ' \2M J M{r/2M-1) 'J ^ ' 

From this equation wc obtain, 

1 + ( - in - l) + M(,^;2M-1) ^^-) = 0- 

Here F = is also the solution of the dynamical horizon. The dynamical horizon 
radius r£) is given by solving (17) as 

re = 2M + 2Me-''/^^. (19) 

Note that this dynamical horizon radius is outside the r = 2M, that is othor 
solution. 
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2.3 Black hole evaporation problem with approximation 

At first, we should derive the energy- momentum tensor T^^ for the integration 
of the dynamical horizon equation. For this end we derive it from the given 
Vaidya matter. For G = 1, F ^l - the non- vanishing components of 

the energy-momentum tensor becomes 

= ^(FA//,-t-M,„) (20) 

Tir' = -^M^ra (21) 

Tr*r' = 0. (22) 

Here we have made the coordinate transformation from r to r* . Writing Tu in 
terms of T^y and T^r* given by (l^ipTj) with 1°- — v"-, we see 

Til -\- Tyj-* 

^ ' ~FMr - M„ - aMr) 



47rr^ 

= ^-M^. (23) 

With being the unit vector in the direction of t", we obtain 

For the dynamical horizon integration we get 

/ ^nrlTi^dro = - (1 + e-^'/''')dM, (25) 
where we have used 

p-v/2M 

F-YT^^^ ('^) 

and the fact 

changing the integration variable from to M. In the above calculation, we 
treat Af ,„ and F^^ with ro fixed, because these functions are used only in the 
integration. Inserting equation (I25p to the dynamical horizon equation ([2]), we 
obtain 



i(2M + 2Me-"/2M) 



Ma 1.M2 g 
Ml -I Ml 2 



(1 + e-'''^^')dM. (28) 
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Taking the limit M2 -> Mi = Af, we obtain 



_^(l + e-"/2A/)^^g-./2M^0^ (29) 
Z Zlvl 

This equation is the dynamical horizon equation in the case that only the Vaidya 
matter is present. There is no solution of this equation, except the trivial one 
(i^ = or r = 2M), so 

ro = 2M{v). (30) 

Here Miv) is the arbitrary function only of the v, which represent the Vaidya 
black hole spacetime. 

Next, we take into account the Hawking radiation. To solve this problem, we 
use two ideas that is to use the dynamical horizon equation, and to use the 
Vaidya metric. The reason to use the dynamical horizon equation comes from 
the fact that we need only information of matter near horizon, without solving 
the full Einstein equation with back reaction being the fourth order differential 
equations, for a massless scalar field. For the matter on the dynamical horizon, 
we use the result of Candelas [5], which assumes that spacetime is almost static 
and is valid near the horizon, r ~ 2M . 

Tti = —Ttt 

1 f°° duoj^ 



2cM47r2(l - 2M/r)' 
where we have used a well known result. 



27r2(l - 2M/r) J„ e^^^^^" - 1 

^ (31) 



1 15a4 ' 



(32) 



and where c = 61440. This matter energy is negative near the event horizon. 
In the dynamical horizon equation, if black hole absorbs negative energy, black 
hole radius decreases. This is one of the motivations to use the negative energy 
tensor. Next we replace length of t to unit length, because in the dynamical 
horizon equation t is used, so 

(33) 

and therefore, the energy tensor becomes 

" 2M4c7r2(l-2Af/r)2' ^^^^ 
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Calculating the integration on the right hand side of ^ for this matter, 



''^ :drn 



M4(l - 2M/rD) 

-bj j^^^, ^dM 

= b j"^ !il±£Z^!^e-"/^^ [2{1 + e-^/^M) ^ ^e-"/2M^ ^^^^ (35) 

Here we insert the expression for r^j (|19p in the first line, and the expression for 
dro/dM = 2(1 + e^''/^*^) + ^e^^/^*^ is used. Here 5 is a constant calculated 
in ^ 

^ (36) 



307207r 



If we also take account of the contribution of the Vaidya matter, and inserting 
this into the integration to the dynamical horizon equation ([2]), we obtain 



2 

M2 r,2 



Ma 
Ml 



Ml 



i(2M + 2Me-^/2*0 

t^}:±I^l!^^e--/M ( 2(1 + e-V2M) iLg-./2M\^^ 
M2 \ ^ ' M J 

(1 + e-"/2A^)dM. (37) 



Ma r 
Ml ^ 



Taking the limit M2 ~^ Mi — M, we finally get 



2^^ ^^2M 



M2 

or 



r , Le-/^' I 2(1 + e-"/^^'^ ) + —e-^r^^' ] , (38) 



- _|(l + e-/2M) + _2_e-/2Ml^^l + ^ )+2m' / ^^^^ 

This is the main result of the previous work that describes how the mass of 
black hole decreases. This equation is the transcendental equation, so usually 
it cannot be solved analytically. However, with the right hand side depending 
only on —v/2AI, we can easily treat Eq. (|39|) numerically. Figure[2]is a graph of 
M as a function of v. By the previous study we found that black hole finally 
vanish at u = and M = 0. The limit of the —v/2M is the solution of the 
following equation as, 

(1 + e'-mi^ ^ ^^-v/2M ^ o_ (40) 
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Figure 2: Numerical calculation of the black hole mass M as a function of v 
from the equation (|39|) 



3 Black hole evaporation without approxima- 
tions 

The result of subsection2-3 use Candela's approximation [5]. However, we 
would like to solve the dynamical horizons equation without approximation. So 
we use stress-energy tensor of massless scalar field. By the similar approach, we 
can obtain following equation as 

i(2Af + 2Me-^^^^') ^ inrl (^cj),ab - |ga6</'x0'") t^^^drc 



M 



+ / -{l + e-'"/^^')dM. 



(41) 



Here, cj) is the massless scalar field. The calculations of the integrand of the first 
term becomes as, 



(42) 
(43) 



Then the first term becomes as 



V ^ / "'0 



27rr^a ^<j),vvdrD. (44) 



Then the dynamical horizons equation 1411 becomes as, 

f-ro i-M 



1 Z"'" 
-(2M + 2Me'''/2A/) ^ 2 

2 Jo 



drD+ / 'Ul + e-'"'''^")dM. (45) 







Using the equation ([27l) and by derivative of M, we can finally obtain following 
equation as, 



(1 + e-) 



f(l + e^) - xe^ 



7r(l + e^)2e-32^ 2(1 + e^) - 2x6" 



(46) 
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Here, x — —v/2M. This equation is main result of our works in this paper. 
This equation is the second order ordinal differential equation with non-linear 
term. This equation can not be solved by analytical method, but it can be 
solved numerical method. The equation is the simple form of the field equation. 

When black hole vanishes, the approximation of the subsection l2.3l is correct, 
then the denominator of the equation becomes 0. And the equation (|46p 
becomes diverge. It means the derivative of the field diverges, and it means 
Schwarzschild spacetime becomes Euclidean spacetime. 

If the X become infinity, the equation becomes simple as 

= (47) 

Next we consider when black hole vanish. Using the following relation 

l + e^^xe^, (48) 
we can simplify the equation (|46p as 

Because of numerical calculation, we solve 

1 + - xe^ = 0, (50) 
then we obtained the x = 1.27846. So the equation (H^ becomes as. 

Here, ai = 4.59108,02 = 46.311. More simplicity we put the right hand side 
of Eq. ((5T|) as constant C, then we can solve this equation. And the solution 
becomes 

(j){x) 2.29554 log[cosh(\/C x 0.308035a;)] . (52) 

For the large C and small x this function behaves like proportional function. 
And when black hole banishes, the function (j){x) becomes 0. 



4 Discussions and conclusions 

We can say black hole finally vanishes by the Hawking effect by our study as 
we first prove [1] . The reason why we enter freedom of the field is that we would 
like to check generalized second low of the black hole thermodynamics. The main 
result of present paper is Eq.ipS]). By this result, numerical calculation of the 
field equation can be simplified. By our study we can show the field closes zero 
in the limit of the black hole evaporation. So we can say that the generalized 
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second low of the black hole thermodynamics is collapsed in this semi-classical 
approach. 

The main problem of our work is at the point that our approach is semi- 
classical one, i.e. the gravity is classical and the field is the quantum one. So 
we should quantize gravitation. And our work may be enlarged to the charged 
black holes and the Carr black hole. 

When the black hole vanishes, the field does not continuous. So we can say 
Schwarzcild universe changes Euclidean universe. 



acknowledgments 

We would like to acknowledge A.Hosoya, M.Shiino and S.A.Hayward for com- 
ments and discussions. 



References 

S.W.Hawking Commun. Math. Phys. 43 199 (1975) 



L.H.Ford gr-qc/9707062 (1997) 



N.D.Birrell and P.C.W.Davies Quantum fields in curved space Cambridge 
University Press (1982) 

S.Sawayama Phys. Rev. D 73 064024 (2006) 
P.Candelas Phys. Rev. D 21 2185 (1980) 
A.Ashtekar, B.Krishnan Phys. Rev. Lett 89 261101 (2002) 
A.Ashtekar, B.Krishnan |gr-qc/0407042| (2005) 
D.N.Page Phys. Rev. D 13 198 (1976) 



C.M.Harris and P.Kanti JHEP 0310 014 hep-ph/0309054 (2003) 



J.G.Russo hep-th/0501132 



T.Tanaka Prog.Theor.Phys.Suppl 148 307 (2003) 

P.Anninos et. al. Phys. Rev. Lett. 74 630 (1995) 

W.A.Hiscock Phys. Rev. D 23 2813 (1981) 

P.Hajicek Phys. Rev. D 36 1065 (1987) 

E.Sorkin and T.Piran Phys. Rev. D 63 124024 (2001) 

R.S.Hamade and J.M.Stewart Class. Quantum Grav. 13 497 (1996) 

I.Brevik and G.Haines Phys. Rev D 67 023508 (2003) 

S.A.Hayward gr-qc/0506126 



11 



